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TRANSIENT ROTORDYNAMIC ANALYSIS 
FOR THE SPACE-SHUTTLE MAIN ENGINE. 

HIGH-PRESSURE OXYGEN TURBOPUMP* 

DARA W. CHILDS** - 
THE UNIVERSITY OF LOUISVILLE 
LOUISVILLE, KY. 

A simulation study has been conducted to examine the 
transient rotordynamics of the space shuttle main engine (SSME) 
high pressure oxygen turbopump (HPOTP) with the objective of 
identifying', anticipating, and (hopefully) avoiding rotordynamic 

problem areas. Simulations were performed for steady-state 

0 

operations at emergency power levels and for critical speed 
transitions., No problems are indicated in steady-state operation 
of the HPOTP at emergency power, levels. However, the results 
indicate that a rubbing condition will be experienced during 
critical speed transition at shutdown. Specifically, for the 
presently projected rotor deceleration rate and imbalance 
distribution, rubbing at the turbine floating-ring seals is 
predicted. This condition can be corrected to some extent by either 
reducing the imbalance at the HPOTP hot-gas turbine wheels, or by 
a more rapid deceleration of the rotor through its critical speed. 

This study was the result of a cooperative effort by the 
author, NASA, and Rocketdyne Division of Rockwell International 
personnel. The rubbing condition indicated for the HPOTP is 
presently under investigation by both NASA and Rocketdyne personnel. 
The HPOTP test program has been designed to determine whether or 
not this rubbing condition will arise, and whether such rubbing 






will compromise the operational performance of the HPOTP. 
Specifically, the test and development program will ascertain- 
the operational integrity of the floating-ring seal. 

This study demonstrates the utility of a transient flexible- 
rotor simulation model as a tool in the development of high-speed 
flexible rotating equipment. Futhermore, the availability of 
the model developed in this study will provide a significant 
analytical reserve for dealing with problems which the SSME may 
see fit to provide once the test program is underway. For 
example, the model would be well suited for the examination of 
alternative hardware "fixes” if an unanticipated unstable whirling 
mode should appear during the development program. 

*This work was performed while the author was a fellow in 
the 1973 NASA-ASEE summer faculty fellowship program at Marshal 

• I 

Space Flight Center. The NASA counter part for this study was 
Mr. Loren Gross. 

**Associate Professor of Mechanical Engineering and 


Engineering Management. 


INTRODUCTION 


A series arrangement of a low-pressure turbopump followed 
by a high-pressure turbopump is employed on both the lox and fuel 
flow paths of the SSME (space shuttle main engine). This paper 
summarizes the significant results of a simulation study (1^)* of 
the transient dynamics of the SSME (space shuttle main engine) 
high pressure oxygen (HPOTP) and fuel (HPFTP) turbopumps, which 
had the primary objective of identifying, anticipating and 
(hopefully) avoiding rotordynamic problems areas. The simulation 
study (_1) consists of a parallel investigation of both the HPOTP 
and the HPFTP. However, since' no problems were revealed for the 
HPFTP, and a potentially serious operational problem is predicted 
for the HPOTP, this paper is restricted to the pertinent results 
for the HPOTP. 

The simulation model used in this study was developed by 
D. Childs (2^) to predict the transient motion of flexible rotating 
equipment. An outline of the structure of this model is provided 
in Appendix A. The model's transient output includes bearing 
displacements and reactions, rotor whirl velocity, rotor spin 
velocity and acceleration, and rotor displacements at. selected 
locations along the shaft. Of these variables, bearing reactions 
and rotor displacement are the most significant in anticipating 
problems which may arise during the operation of the turbopumps. 
Excessive bearing loads result in shortened bearing lifetimes, 
and excessive deflections cause seal failures; and can potentially 
lead to catastrophic failure m.odes . 


^Denotes listing in reference section. 


The SSME turbopumps are presently under development by 
Rocketdyne Division of Rockwell International, Canoga Park, 
California. The basic design of the engines is fixed with the 
development of detailed drawings in progress; however, component 
testing of turbopump hardware will not begin until approximately 
June 1974. Rocketdyne 's analysis which is pertinent to a 
rotordynamic investigation is largely complete, and the steps 
involved in this analysis j together with the requisite results, 
are listed below: 

a. Structural dynamic analysis of the’ SSME rotating assemblies. 
This analysis yields a structural dynamic model for the rotating 
assembly, which is customarily used to calculate rotor Critical 
speed locations and whirling modes. As an aid to the present 
investigation, Rocketdyne personnel also calculated unsupported 

or "free-free” natural frequencies and bending modes for the rotor. 

b. Calculation of load-deflection properties for the SSME 
ball bearings. A canputer program which is based on the analyses 
of A. B. Jones (_3 ) , (4) is used to calculate the radial stiffness 
of ball bearings as a function of the axial and radial load 
carried by the bearing and the running speed of the shaft. 

c. The magnitude of the hydrodynamic side loads which act 
on the turbopump rotating assembly have been estimated. 

d. The magnitude and axial distribution of rotor imbalance 
have been estimated. These estimates are based on both the planned 
balance procedure and past experiences with similar turbopumps. 
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The data cited are sufficient to define a transient rotor- 
dynamic model. The assistance of Rocketdyne personnel in pro- 
viding this information is greatly appreciated. In particular, 
the efforts of Mr. Bernie Rowan, Mr. Miles Butner, and Mr. Bill 
Brooks have been of considerable value. 

Operational testing of the SSME turbopumps will begin in 
October 1974. The instrumentation which is planned for the 
test phase will include measurements of transient radial 
motion of the rotating assembly at five axial locations. The 
simulation study (1^) provides the transient motion for both 
turbopumps at these locations, and should assist in the under- 
standing of the transient motion data when it becomes available 
Conversely, the availability of high-quality data (l.e., 
calibrated, high-resolution, transient data taken along the 
length of the rotating assembly) will provide a unique 
opportunity to remove many of the uncertainties associated with 
transient rotordynamic modeling. In particular, this data, in 
conjunction with the. present simulation models, will enable a 
much better indentification of the damping provided by ball 
bearings, transient and steady-state hydrodynamic forces, rotor 
imbalance distribution, etc. 

The availability of the two models developed in (1.) will 
provide a significant analytical reserve for dealing with 
problems which the SSME turbopumps may see fit to provide once 
the test program is underway. For example, the models would 


( 3 ) 


be well suited for the examination of alternative hardware "fixes" 
if an unanticipated unstable whirling mode should appear during 
the development program. 
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SPACE SHUTTLE MAIN ENGINE HIGH PRESSURE OXYGEN TURBOPUMP 

. • 

General Information . 

The SSME HPOTP rotating assembly is illustrated in figure 1. 

At normal power levels, t*he main stage of this turbopump is to be 
supplied lox by a low pressure turbopump at a flowrate of approxi- 
mately 1047 lb. /sec. at a pressure of approximately 360 psia. The 
discharge pressure of this stage is to be approximately 4635 psia, 
and the bulk of the discharge flow goes directly to the thrust 
chamber; however, approximately 100 lb. /sec. is diverted to the 
boost stage with an inlet pressure of a'pproximate ly 4390 psia. 

The discharge flowrate from the boost stage (7680 psia) feeds the 
preburners which drive the hot gas turbines for both high-pressure 
turbopumps. The gas flowrate which drives the turbine approaches 
it from the rear, i.e., it moves axially towards the pump stages. 

The turbine develops approximately 22,900 horsepower at the nominal 
running speed of 29,250 rpm. Turbine torque at these conditions is 
4060 ft. lbs. Discharge flowrate from the turbine goes directly 
to the main thrust chamber. At emergency power levels, the turbo- 
pump running speed is increased to 31, 160 rpm, while* at minimum 
power level, the speed is reduced to 20,890 rpm. 

Rocketdyne analysis (^) indicates that a hydrodynamic side- 
load is developed on the rotating assembly of this turbopump due 
to circumferential pressure variations on both the main and boost 
stages. These pressure variations are the result of single-discharge 
designs for both stages, and are assumed here to be proportional 
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to ^ (spin-velocity squared) . They are modeled by 

— 5 ^ 2 . 

(boost) = ^ /.go^/o /i- 

f(main) = ^ ^ / . S'J X /iS . /^- (1) 

# 

Hence, at emergency power levels, (p = 31,160 rpm (2990 rad. /sec.) 
the side loads for the boost and main stages are 161 and 1370 lbs., 
respectively. The direction of the loads is fixed relative to the 
turbopump housing. The two loads are oppositely directed. 

The bearings in the HPOTP are designed to react radial loads 
only. A balance-piston arrangement is provided at the discharge 
of the mainstage to react axial thrust loadings. The ball bearing 
used for the forward^duplex-bearing set of the HPOTP is the same 
as that used for both bearing sets in the HPFTP. The forward and 
rear bearings preload for the HPOTP are 350 and 550 lbs., respectively. 
This preload is maintained by loading the outboard and inboard 
bearings against one another via a Belleville spring. 

As would be expected, the radial load-deflection characteristics 
of a ball bearing depends on both the axial and' radial load carried 
across the bearing and the bearing's speed of .rotation. For the 
magnitude of radial bearing loads developed in this turbopump, the 
relationship between radial loads and deflections is linear, ^nd 
is appropriately modeled by a bearing stiffness coefficient. The 
dependence of this bearing stiffness coefficient on running speed 
and axial preload is illustrated in figures 2 and 3. The data of 
figure 2 are for the rear bearings of the HPOTP, while the data 
for figure 3 are for both the forward bearings of the HPOTP and 
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both HPFTP bearings . These data were generated by Mr. M. F. 

. • 

Butner, Rocketdyne Division of Rockwell International, Canoga 
Park, California, via a computerized version of A. B. Jones analyses 
O ) , (A). Since a constant axial preload is a design objective 
in the bearings, radial s*tiffness is only dependent on shaft spin 
Speed. For the data presented, the following representation is 
employed 

3 ^/ ^ /OS'. 2 /j 6. 

* "". 5 / 2 ^ /jb. 

sr. 3/2 X JO^ SCC?/^ . 

kj^ CK(rer ) = f 

3ro - i^'32/ X/0^ 

3r/ = BA Z0 S^C. /c^ . 

^rz - - /esr AA>. f Ah . 

S ^ ^ ' Cf S Z Zo • S'(ZC ’ . ^ . 
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Damping coefficients used in this study (arbitrarily) to account 
for dissipation at the forward an,d rear bearings are 

• 

- , ‘ (4) 

SCC' / • 

% 

A principal objective of the bearing design in the HPOTP is 
critical-speed- location control. The first two -free-free bending 
modes (zero-bearing stiffnesses) for the rotating assembly of this 
turbopump are located at approximately 25,570 and 58,170 rpm. 

Because of the comparatively* stiff bearing design (and gyroscopic 
stiffening), the calculated critical speeds of the turbopump lie 
at approximately 13,000 and 40,000 rpm. Hence, the operating speed 
range of the turbopump (20,890 to 31,160 rpm.) is well removed from 
the ro.tor's critical speeds. 

The balance procedure presently planned for the HPOTP consists 
of the following steps: 

a. The rotating assembly is balanced at low speed in.'.a balance 
assembly. 

b. The rotating assembly is then partially disassembled prior 
to reassembly in the turbopump. 

On the basis of past experience with similar turbopumps, Rocketdyne 
personnel* suggest that reasonable values for residual turbopurap 
imbalance would be the following: 

a. 1 gm. in. at the boost stage impeller, 

b. 5 gm. in. at the main impeller, and 

c. 5 gm. in. at the turbines. 


*Personal communication with B. Rowan, 6/2/73. 
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The orientation of the rotor imbalances relative to either one 

. • 

another or the rotor is not predictable. 

At the present time, Rocketdyne plains to measure the transient 
motion of the rotor with five radially mounted and two axially 
mounted Bently proximity gauges. The planned radial Bently 
locations are as follows: 

a. Immediately behind the boost stage impeller, and 
approximately 2.6 in. from the front of the rotating assembly, 

b. Immediately behind the inboard front bearing, and 
approximately 5.4 in. from the front of the rotating assembly, 

c. At the forward inlet of the main stage impeller, and 
approximagely 7.6 in. from the front of the rotating assembly, 

d. Immediately forward of the inboard rear bearing, and 
approximately 13.9 in. from the front of the rotating assembly,, and 

e. Forward of the inboard hot -gas turbine, and approximately 
20.13 in. from the front of the rotating assembly. 

The axial Bently's are to be located at the fear of the boost and 
mains tage impellers. 

HPOTP Structural Dynamic Properties 

The structural dynamic analysis for a turbopump is initiated 
by modeling its rotating assembly as a collection of rigid bodies 
connected by a massless elastic structure. Rocketdyne 's structural 
dynamic model of the HPOTP employs 13 rigid bodies connected by 
either uniform or discontinuous beam structures. Each rigid body 
has both a displacement and a rotational degree of freedom. The 
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inertia properties of the rigid bodies, together with their 
location relative to the mass center of the rotating assembly., 
are provided in Table 1. The entries are the mass, 
diametral, . and polar moments of inertia of the ith component 
rigid body . The entries in Table 1 are the coordinates 

of the component, rigid bodies in the coordinate system. The 

origin of the system coincides with the rotating assembly 

mass center. The component defines the axial position of the 
rigid bodies with positive ^ being directed towards the turbine. 
The and components define the eccentricity (imbalance) 

of component rigid bodies. The magnitude and axial location of 
the imbalance distribution defined in Table A.l is consistent 
with Rocketdyne estimates; however, the orientation of the im- 
balance with respect to the system is arbitrary. The 

rigid body properties of the HPOTP follow 

^9,9 (,3 /i . s'ecf. V 
^ ~ ^ ^^2. /6. sec ^ 

^ -Sr.J04 K /A 


( 5 ) 


/b. sec. 




-J 


-J . 
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where m is the rotor mass and J and J are the rigid body diametral 
and polar moments of inertia. Further, the components 3^ j 3^ 
define rigid body imbalance and ^ are the rigid body 

products of inertia which result from the imbalance distribution 
of Table 1. 

The results of eigenanalysis performed by Rocketdyne which are 
pertinent to this investigation are the lowest two (non-zero) 
free-free eigenvalues and their associated eigenvectors. The two 
lowest free-free natural frequencies follow 


^ 2jC.77.3 - 

- C/Oe>/ .9 jre£f. ^ 3 , 2.02 


An eigenvector for the rotor includes both the displacement and 
rotational degrees of freedom, for each component body. These 
eigenvectors are presented in partitioned form in Table 2 where r/leJ 
contains the displacements and the rotational degrees of 

freedom. The eigenvectors are normalized to satisfy the following 
orthogonality condition 


where and are diagonal matrices whose elements are and 

Vi . The superscript T in Eq. (7) defines the matrix transpose 
operation , and iTTl is the unity matrix. Stations 2, 6, 8, 11, 
and 17 of Table 2 correspond to axial positions for which Bently 
proximity gauges are planned. . These station locations were not 
included in Rocketdyne ’s original structural dynamic model. The 
entries in Table 2 for these station locations were obtained via 
linear interpolation between adjacent stations. 
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1. 

-12.75 

8.67 

10.61 

16.92 

0.659 

0.0 

2. 

-11.70 

6.36 

4.672 

8.284 

0.0 

0.0 

3.- 

-10.69 

2.24 

0.833 

1.314 

0.0 

0.0 

4. 

-9.82 

1.79 

0.618 

1.055 

0.0 

0.0 

5. 

-8.19 

1 Ml 

5.976 

6.620 

0.0 

0.0 

6. 

-5.02 

62.1 

• 228.2 

257.7 

0.0 

0.460 

7. 

-1.42 

9.34 

8.967 

8.560 

0.0 

0.0 

8. 

.126 

3.55 

1.927 

1.927 

0.0 

0.0 

9. 

1.25 

5.74 

3 . 666 

3.665 

0.0 

0.0 

10. 

2.31 

6 . 64 

7.213 

7.212 

o.on 

0.0 

11. 

3.75 

7.81 

8.079 

8,079 

0.0 

0.0 

12. 

5.A8 

4.69 , 

30.99 

3.099 

0.0 

0.0 

13. 

7.63 

70.2 

508.8 

508.8 

0.407 

0.0 



TABLE 

1. Geometric and : 

inertial data for the SSME HPOTP 


3*- 

mer^r}yt 



1. 

-12.75 

5.112 

,3.965 

1.025 

1.495 

2. 

-12.06 

4.401 

2.915 

1.023 

1.486 

3. 

-11.70 

4.032 

2.369 

1.022 

1.480 

4. 

-10.69 

- 2.983 

.8182 

.9998 

1.384 

5. 

-9.82 

2 . 106 

-.4005 

.9548 

1.199 

6. 

-9.25 

1.595 

-.9250 

.8628 

.8493 

7. 

-8.19 

/6602 

-1.886 

.6943 

.2077 

8. 

-7.00 

-.0528 

-1.877 

.6401 

.0384 

9. 

-5.02 

-1.221 

-1.861 

.5514 

-.2391 

10. 

-1.42 

-2.771 

.0513 

.0381 

-.5219 

11. 

-.749 

-2.698 

.4687 

-.0598 

. -.4980 

12. 

.126 

-2.604 

1.009 

-.1865 

■ -.4671 

13. 

11.25 

-2.340 

1.573 

-.2738 

-.3709 

14. 

2.31 

-1.998 

1.972 

-.3523 

-.2337 

15. 

3.75 

-1.397 

2.236 

-14534 

.0223 

16. 

5.48 

-.4504 

1.921 

-.5891 

.4973 

17. 

5.50 

-.4304 

1.903 

-•6090 

.6428 

18. 

7.63 

.9311 

.377 

-.6614 

.9092 


TABLE 

2 . HPOTP 

rotating 

assembly eigenvec,t;ors 
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Simulation Procedure and Results 


The output capability of the simulation model is 

summarized below. 

1. Bearing motion. The transient motion of the shaft at 
bearing locations is generated. Bearing motion output 
includes both the motion of the shaft in two orthogonal 
planes as a function of time and X-Y plots of shaft orbital 
motion, i.e., the output corresponds to the shaft motion 
which would be observed from either a side or axial view 

of the rotor. Solutions are generated for the forward- 
inboard and rear-outboard bearings. 

2. Modal coordinates. The transient solution of modal 
coordinates associated with the two free-free rotor bending 
modes are generated. These solutions are of inferest in 
verifying the correctness of the model's operation, but 
are not significant from a turbopump operational viewpoint. 

3. Bearing reactions. The resultant reactions for forward- 
inboard and rear-outboard reactions are generated as a 
function of time. 

4. Spin velocity. The rotor spin velocity in rpm. is generated 
as a function of time. 

2 

5. Spin acceleration. The rotor spin acceleration in rad. /sec. 
is generated as a function of time. 

6. ' Simulated Bently output. Motion of the rotor is generated 

for axial positions along the shaft corresponding to planned 
Bently proximity gauge locations. Radial motion in two 
orthogonal planes is defined with both time' and X-Y solutions 


0.3) 


7. Maximum shaft deflections. Plots of radial shaft-deflection 
magnitudes as a function of rotor axial position are generated 
periodically. The following procedure is followed in obtaining 
these plots. At each *time tj^^ in the simulation, a vector is 
generated whose components are the shaft-deflection magnitudes 
Wj(t^) at j rotor axial stations. The maximum of the j rotor 
deflection magnitudes Wj (t , denoted W(t^) , is calculated. 

A search is made over 100 succeeding values of tj_ for the 
maximum value of W(t£), denoted . The plots of 

maximum shaft deflection magnitudes provided correspond to 

•IL * 

the Wj (t'j^) vectors. In other words, a plotted shaft- 
def lection-magnitude vector includes the maximum rotor 
deflection magnitude over the time of observation. 

The validity of a transient simulation model is primarily 
established by the degree to which the model correctly identifies 
the location and form of rotor critical speeds. As noted previously 
the critical speed locations calculated by Rocketdyne for the 
HPOTP rotating assembly are approximately 13,000 and 40,000 rpm. 

The transient model correctly identifies both the critical speed 
locations and the deflection mode shape of the rotor during 
transition through the critical speed. 

From a rotord 3 mamic viewpoint, the most severe operating 
conditions for the HPOTP arise during transition through the 
first critical speed and during steady-state operation at 
emergency power levels. Considering the, critical speed transition 
first, the amplitude of transient rotor deflections largely 
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depends on the rate at which the rotor is accelerated ( or de- 
celerated) through the critical speed and the damping which is 
present. The proposed system operation of the SSME is such 
that the critical-speed transition during startup is considerably 
more rapid than at shutdown; hence, the shutdown transition was 
examined. The rate of deceleration of the rotor in the vicinity 
of the first critical speed (^) during shutdown is approximately 
940 rad. /sec. while the comparable value during startup is 

O 

3660 rad. /sec. 

Damping is provided in the rotor model by modal damping in 
the rotor arid by concentrated damping at the bearings. The 
modal damping factors used in this study were two percent of 
critical. This value is felt to be conservative. Damping at 
the bearings is accounted for by linear damping coefficients 
which would correspond to the physical model of a dashpot 
attached between the turbopump housing and the rotating assembly. 
There is a great deal of uncertainty involved in the amount of 
damping provided by the ball bearings and the fluid surrounding 
the turbopump rotating assembly. The damping factors used in 
this study yield a rapid attentuafiOn of initial starting 
transients, and are felt to reasonably account for external 
rotor damping . 

Figure 4 illustrates some of the results of a transition 

through the 13,000 rpm. critical speed. The rotor is decelerated 

from 13,380 to 12,500 rpm. at a nominal deceleration rate of 
2 

940 rad. /sec. . The results of the separate frames of figure 4 
are discussed below. 
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Figures 4(a), (b), and (c) illustrate the motion of the ^ 

shaft at the turbine floating-ring seal location. Figures (a), 
and (b) illustrate the shaft-deflection time history in the 
orthogonal X-Y and Y-Z planes where the Z axis defines the 
nominal geometric axis of the turbopump, i.e., these solutions 
correspond to the motion which would be observed by looking 
in at the shaft radially. By contrast, the solution presented 
in figure (c) corresponds to the motion which would be observed 

tf. 

from an axial position. The steady-state biased displacement 
of the rotor displayed in figures 4(a) and (c) results from 
the hydrodynamic side loads, which are applied to the turbopump 
in the X-Z plane. 

Figures 4(d) through (h) illustrate the maximum shaft 
magnitudes experienced during initial-speed transition. A 
better appreciation of this figure is obtained by noting that 
the magnitude (absolute value) of deflection is plotted. The * 
markings on the plot denotes either a station location in Rocketdyne's 
structural dynamic model or a planned Bently location. The forward 
bearings on the turbopump are located at. -10.7 and -9.8 inches, 
the rear bearings are located at .12 and 1.24 inches, and the 
extreme right-hand station corresponds to a location between 
the turbine wheels. The figure basically illustrates that 
large deflections are the result.^:of the overhung turbine wheel 
design. The deflections in the neighborhood of 5.5 in. are 
sufficient to cause rubbing on the turbine floating ring seals. 

The clearance at this location is .005 in. and the predicted 
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deflections are .010 in. Time solutions for the deflections at 
this location were provided in figures 4(a), (b), and (c). 

The operating conditions and physical -parameters which were 
used in the simulation model to obtain a rubbing condition at the 
turbine-floating -ring seals are not worst case. Specifically, 
the imbalance distribution, the damping and the time interval 
required for critical-speed transition could all be worse than 
those used here with a consquent increase in deflections. The 
imbalance distribution of Table 1 was not chosen to maximize the 
rotor deflections. ^ For example, larger deflection amplitudes 
could be obtained from the following distribution 


with all other , a;r zero. The two percent of critical 

>J 

damping factors used to model structural damping of the rotating 
assembly could conceivably be half percent or less. A simulation 
run was made in which a 100-lb. step increase in side loading was 
imposed on the rotating assembly (at zero running speed). The 
resultant damping factor was approximately 3.5 percent of critical. 
The same type of run was made at a running speed slightly above 
the critical speed location, and the motion due to the disturbance 
was completely absent from the plots after 4 cycles demonstrating 
that the model is not lightly damped. The deceleration rate 
through the critical speed used here was 940 rad. /sec. which 
is the predicted (6) rate between 14,000 and 13,000 rpm. However, 
the predicted deceleration rate for less than 13,000 rpm. is 
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approximately 520 r^d./sec. Hence, the deflections predicted 
in this stjudy are considered to be reasonable, and a possibility 
exists for larger and potentially more damaging amplitudes. 

The transient response of the modal coordinates 
are illustrated in figures 4(i) and (j). The peak magnitudes 
of these functions coincide (in time) with the shaft motion of 
figures 4(a) and (b) . The amplitude of is seen to be an 
order of magnitude larger than that of elastic 

d . 

deflections of the rotor are defined in the rotor-fixed x, y, z 
coordinate system as a linear summation of the modal coordinates. 

The oscilliation in these coordinates at the running speed 
frequency results from the hydrodynamic side loads. The combined 
effects of rotor imbalance and the critical-speed transition 
account for the comparatively slow oscillation in these variables. 

The rotor spin acceleration is illustrated in figure 4(k). 

This solution is included to demonstrate the coupling of the 
transverse oscillations to the spin-axis motion. This feature 
of the simulation model is of no particular value in the present 
(high-spin acceleration) application. However, it is required to 
correctly simulate a "slow", power- limited, cirtical-speed transition. 

Although not illustrated, the peak reactions for forward and 
rear bearings are approximately 220 and 450 lbs., respectively. 

This magnitude of transient bearing loads does not constitute a 
pr ob lem . 

A simulation run was .executed to examine the motion of the 
turbopump at EPL, which corresponds to = 31,140 rpm. The 
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results basically resemble those presented in figure 4. Figure 
5 illustrates a shaft deflection magnitude plot for EPL conditions, 
and demonstrates, that shaft deflections are generally less than 
.002 in. The maximum bearing loads on the forward and rear 
bearings at EPL are approximately 750 and 650 lbs., respectively. 
The results of the HPOTP simulation runs for EPL conditions do 
not indicate any problem areas associated with either excessive 
shaft deflections or excessive bearing loads. 
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CONCLUDING COMMENTS 


The results of this study indicate no problems in operating 
the HPOTP at EPL conditions. However, the results indicate a 
potentially serious problem during critical-speed transition at 
shutdown. Specifically, for the presently projected rotor 
deceleration rate and imbalance distribution, a rubbing condition 
at the turbine floating-ring seals is predicted. The problem can 
be reduced to some extent by reducing the imbalance at the HPOTP 
turbine wheels. The problem would also be reduced if the rotor 
were decelerated through the critical speed more rapidly, or if 
additional damping were provided. , 

This study resulted from the cooperative assistance of both 
NASA and Rocketdyne personnel, and (to the author's knowledge) 
represents the first use of a general purpose, transient, 
flexible -rotor simulation model as a development tool in a 
liquid-rocket-engine program prior to hardware testing. This 
study was conducted to anticipate and avoid hardware failures. 

By contrast, simulation efforts in the past have resulted from 
(sometimes catastrophic) hardware failures, and have been hampered 
by ad hoc models of limited capability. The results of this study 
have indicated a potentially serious operational problem well in 
advance of operational turbopump testing. Time has accordingly 
been purchased for a careful analysis of alternative hardware and 
operational adjustments. Time has also been purchased in the 
design of the turbopump test program to verify the effectiveness 
of such adjustments. 
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APPENDIX A 


THE SIMULATION MODEL - ' 

Kinematics 

A modal simulation model from Ref. [ 2 ] was employed in this 
study. The rigid body kinematics required are illustrated in 
figure A.l. The x, y,z coordinate system of this figure is fixed 
in the elastic rotor, and its origin coincides with the rotor's 
mass center. The z axis coincides with the nominal rotational 
axis of symmetry of the rotor. The vector R locates the mass 
center of the rotor in the inertial X, Y, Z system, and fi is the 
angular velocity of the x, y, z system relative to the X, Y, Z 
system. The orientation of the x, y, z system relative to the 
X, Y, Z system is defined by the Euler angle set , <P> The 

constraints of the bearings on rotor motion justifiy the custo- 
mary' "small" angle assumption for and a^, i.e. , 

= (n-c<{) - n-S({))/ca2= - n-s(j) 

X y X y 

*2 = (A.l) 

i ” diSa, = fi- - 

^ z * 2 z 

where c<j> = cos<() , S(() = sin<j> , etc. 

The operation of differentiation with respect to time in the 
various coordinate systems is defined in terms of the arbitrary 
vector V as follows : 

■O' = y Y 7 “ time rate of change of the vector v 

' ' as seen by an observer fixed in inertial 

space 

V = ^ 1 ^^ y 2 “ time rate of change of the vector v_ 

as seen by an observer fixed in the x, y, 
z system. 


The components of vectors are required in both the X, Y, Z and 
X/ Yr z system, and the following notation is defined in terms 
of the arbitrary vector v 

(v)^; implies components in X, Y, Z system 
• (v)j^; implies components in x, y, z system 
The following orthogonal coordinate transformation relating these 
components is required: 

(v)^ = [A] (V) j , 

where to a comparable, level of accuracy with Eq. (1) 


c<f) 

S(|) . 


+ajS<j) 


-S(j) 

C<t> 

ajscj) 

+aiC<j) 

(A. 2) 

“2 


1 




The rotor's distributed parameter character is accounted 
for by an n-body lumped parameter approximation. The position of 
the ith rigid body relative to the x, y, z coordinate system is 
defined by the vector r^ illustrated in figure A .2 . The vector e^ 
defines the elastic displacement of the origin of the x^, y^ , 
system.. The angular orientation of the body-fixed x^ , ^ 

system relative to the x^, y^, z^ system (and, x, y, z system) is 
defined by the Euler angle set 3^-, analgous in 

their order of rotation to the Euler angle set (j) illustrated 

in figure A. 1(b). Axial and torsional elastic deflections are 
neglected in this analysis; hence, the elastic deflections of the 
rotor are defined in the x-z and y-z planes by the vectors 


(e;)^ = (ei. 

e-^ , 
x' 

••®E> 

- (Sj)’' = 

to 

...3^ ) 


,e-2 , . 

y 

-•®p) 

, (Sj)’’ = 

/ 3j,. 

n (A. 3) 

...3| ) 


where the superscript T denotes the matrix transpose operation. 

The free-free bending modes (eigen-vectors) of the rotor are 
used to define the following orthogonal transformation and to 
introduce modal coordinates 

(e-) = [Ag] (q-) , (e-) = [Ag] (q-) 


(A. 4) 


(e-) = [Ag] (q-) , (e-) =-[Ag](q-) 

The matrices [Ae ] , [Ag] satisfy 

[Ag]'^(m) = 0 , [Ag]T(ma-) + [Ag]'^(J) = 0 (A. 5) 

In addition, they are normalized to satisfy the orthogonality 
relationship 


[Agj'^(m) [Ael + [Ag ] 'T [ j] [Ag ] = [U] (A. 6) 

where [U] is the unity matrix. In Eq. (A. 5), the elements of the 
vectors (m) , (ma-) , and (J) are m^, mia^-, and J^, respectively, 
where 


mj^ = mass of component rigid body i. • 

= diametral moment of inertia of component rigid body i . 
a^ 2 ' nominal axial position of the mass 'center of component 
rigid body i. 

The output variables of primary interest from a rotor simula- 
tion model are the ‘displacements w^ and velocities wi of the 
bearings and of component rigid bodies. The displacement of a 
component rigid body can be defined by the vector equation 
w^ = R - b^ + ai + e^ 

where b^ locates the initial position of the ith component rigid 
body in the X, Y, Z system, and a^ locates its nominal position in 
the x, y, z system. The vector definition for the velocity of 
the rigid body in the X, Y, Z system is 



= R + ei + n X (a^ + e^) 

The component statement of these equations is 
(wi)j = (R)i - (bi)j + [A]T(ai + e^hj 

(wi)l = (R)i + [A]T{(ei)j + t(n)j](ai + e^)^} 

where the notation implies 


(A. 7) 


[(fi)j] = 


0 

-fl- 

n- 


z 

y 

ft- 

0 

-n- 

z 


X 

-Q- 

fi- 

0 

y 

X 



and performs the matrix equivalent of the vector cross product 
operation. From Eq. (A. 4), one obtains the definitions 


(e-) = [Ae] (q-) , (e-) = [Ag] (q-) 


(e-) = [Ael (q-) , (e-) = [Ag] (q-) 


(A. 8) 


Eqs. (A. 7) and (A. 8) can be used to calculate the displacements 
and velocities of each component rigid body. 

This completes the required kinematics for the governing 
equations of motion. The vectors R and n, and the Euler angles 
<* 1 / 012 / used to define the location and orientation of the 

X/ y, z coordinate system, while the modal coordinates (q-) , 

(q-) define the elastic deflections of the component rigid bodies 
relative to the x, y, z system. 


Equations of Motion 

The modal vibration equations for the x-z and y-z planes are 
Cq-) = (F-) + ^{2[U]-lCM] }(q-)-{ [A] + ^2[CK]}{q-) 

+$[CK](q-) + ifi-(Cv) + (GF)- 
(q-) = (F-) - <J{2[U]-[CM] }(q-)-{ [A] + [CK] } (q-) 

+4>’lCK](q-) + ^fi-(Cv) + (GF)- 

X. y y 


(A. 9) 



In the above, [A] is the diagonal matrix of free-free (non-zero) 
eigenvalues. The vectors (F-) and (F-) are generalized modal 


force vectors defined by 


(F-) = 
X 

[Ag]’' 


+ lAei'^(Mj) 



(Ff) = 

lAgl’' 


- [Aj)'^(Mj) 



where 






- = 


4 ' 

... t x) / 

(Mj)’’ = (Mj, m|, . 

. . M*^-) 

y 


< 4 / 

4 ■ 

... f”-) , 
y 

(Mj)'f- (Mi, Mi, . 

. . M^-) 

X 

The vectors 

(CF>- 

and 

(GF)- define 

the modal forcing 

functions 


due to distributed imbalance and asymmetry, and are defined by 


(GF)- = [Z](gf)- , (GF)- = [Z](gf)- 


(A.IO) 


where 


12 


(gf)“ = <!> 

«« 

<!> 

♦ 


(gf)- = -^2 

• * 


and 


[Z] = [(zM (z2) (z3) (z*^)] 

= [ [Ael'^Cma-) , [Agl'^Cma-), (Ag]'^(j2-), [Ag]’^(j2-)] (A.ll) 

Ox ,-rO 


The components of the vectors (ma-) , (ma-) , (J--) , (Jr~) are 

" y -- - X2 




ia^-, nij^a^-, J", J", respectively, where 


a ^-7 a^- = eccentricity of component rigid body i. 

X y 

J", J~ = product of inertias of component rigid body i in 
the body-fixed xi, y^, z± coordinate system. 


Returning to Eq. (A. 7), 

[CM] = [Agl'^EJ] [Ag] 

[CK] = [A3]^[j-J] [Ag] - [Ael'^[m] [Ael = [CM] - [U] 
(c^) = [Ag]'’‘’(j) , (j)T = (jl, j2^ jnj 


(A. 12) 


where [J] and [J-J] are diagonal matrices whose entries are Jj_ 
Xpolar moment of inertia of component rigid body i) and 
respectively. 

The defining equations of motion for the "rigid-body" variables 
R, n are 


(F^)j = m(R)i (A. 13) 

' (M^) ■; = [J];](«)j+ [(fi) j] [J°j] (A. 14) 

l(cv)'^{(qx) + ^(qy)> 

+ l:(c^)'^{(q-) - 

where and are the resultant external forces and moments 
acting on the body, and 

P- = 2| { (zM'^+(z3)^}(q-) + 2^{ (z2)T+(z'+)^} (q-) 

+JJ-(Cv)'^Hq-) + (q~) } + n-(Cv)'^{ (q-)-'J (q-) i ' (A. 15) 


Returning to Eqs. (A. 13) and (A. 14), m is the total mass of the 


rotor, [J j] is the aggregate "rigid body" niomeut of inertia, for the 
rotor in the ;x, y, z system, and 
,o 


iJjJ = [J°j] + 


(Cv)'^(q-) 

(Cv)'^(qy) 


(o^)’'(qj) (Cv)’'(qj> 2M" 

5 J 22 = + (z^)''^} (q-) + { (z2) + (z**) ^ (q-) 


(A. 16) 
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